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Abstract Quantum Optomechanical system serves as an interface for coupling be¬ 
tween photons and phonons due to mechanical oscillations. We use the Heisenberg- 
Langevin approach under Markovian white noise approximation to study a quadrat¬ 
ically coupled optomechanical system which contains a thin dielectric membrane 
quadratically coupled to the cavity field. A decorrelation method is employed to 
solve for a larger number of coupled equations. Transient mean number of cavity 
photons and phonons that provide dynamical behaviour are computed under differ¬ 
ent coupling regime. We have also obtained the two-boson second-order correlation 
functions for the cavity field, membrane oscillator and their cross correlations that 
provide nonclassical properties governed by quadratic optomechanical system. 


1 INTRODUCTION 


Recent years have seen significant experimental progress in realizing deterministic 
interactions between single photons, which has profound importance for future op¬ 
tical technologies. Most novel experiments have been explored in cavity quantum 
electrodynamics (cQED) [T]|2], where photons inherent the saturation of a single 
two-level atom due to strong interactions between the atom and the cavity field. 
A single atom-cavity system described by well-known Jaynes-Cummings model [3] 
has been used as an important test bed for implementation of quantum informa¬ 
tion algorithms as well as construction of a quantum network with the aim for 
quantum computation. The strong-coupling regime of the cQED is reached when 
the coupling strength of the atom with the cavity mode dominates over the de¬ 
coherence processes, due to spontaneous emission from the excited atomic level 
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to ground state and the leakage of photons in the cavity mode. Based on Fabry- 
Perot interferometry [4], semiconductor microcavities have been developed where 
excitons act as quantum systems [5]. Alternative approaches have been explored 
based on slow-light enhanced Kerr nonlinearites [6] , single dye-molecules [7] , strong 
photon interactions mediated by Rydberg atoms [8] , atoms coupled to hollow-core 
fibers [9] and nanofiber traps m- 

Optomechanical system provides a more promising approach towards realiz¬ 
ing strong photon interactions m- In comparison to atomic systems, for example 
macroscopic mechanical resonators are relatively easy to control and to be inte¬ 
grated with other systems. In these systems, an optical cavity, with a movable 
end mirror msM or with a micromechanical membrane is subjected to mechan¬ 
ical effect caused by light through radiation pressure M- So, Cavity Quantum 
Optomechanics has emerged as an interesting area for studying quantum features 
at the mesoscale, where it is possible to control the quantum state of mechani¬ 
cal oscillators by their coupling to light field m- Recent advances in this area 
includes the realization of quantum-coherent coupling of a mechanical oscillator 
with an optical cavity m. where the coupling rate exceeds both the cavity and 
mechanical motion decoherence rate, laser cooling of a nanomechanical oscillator 
to its ground state HU- 

Experimental and theoretical proposals aiming to study quantum aspects of 
interactions between the optical cavities and mechanical objects have focused on 
cavities in which one of the cavity’s mirror is free to move for example, in response 
to radiation pressure exerted by light in the cavity. This nearly includes all op¬ 
tomechanical systems described in the literature, including cavities with ’folded’ 
geometries, cavities in which multiple mirrors are free to move m, and whisper¬ 
ing gallery mode resonators m in which light is confined to a waveguide. All 
these methods have two important considerations. First, cavity’s detuning is pro¬ 
portional to the displacement of a mechanical degree of freedom, which is mirror 
displacement or waveguide elongation. Second, a single device must provide both 
optical confinement and mechanical effects. So, achieving good mechanical and 
optical properties both simultaneously has been the main aim in quantum op¬ 
tomechanical systems because high-finesse mirrors are not easily integrated into 
micromachined devices. Experimental work done so far have achieved sufficient op¬ 
tomechanical coupling to laser-cooled mechanical devices Esun] but have not yet 
reached to quantum regime due to technical problems mentioned above. A more 
fundamental challenge is to read out the mechanical element’s energy eigenstate. 
Displacement measurements cannot determine an oscillator’s energy eigenstate 
and measurements coupling to quantities other than displacement have been dif¬ 
ficult to realize in practice [221123] . 

In this work, we theoretically study an optomechanical system realized in sem¬ 
inal experimental works p4ll^ and has potential to resolve both of these above- 
mentioned challenges. We study a system which contains a thin dielectric mem¬ 
brane with quadratic response to the cavity fields. Coupled Heisenberg-Langevin 
equations are obtained under Markovian white noise approximation, and solutions 
up to second-order correlation operators have been developed in Section II. In the 
same section, we also discuss the transient dynamics of the system Hamiltonian 
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in various scenario. Section III discusses the results with the nonclassical photon 
statistics of the cavity mode as well as mechanical oscillations mode including 
cross-correlation between them. We conclude our results in Section IV. 


2 Heisenberg-Langevin Approach For Optomechanics 


We consider a quadratically coupled optomechanical resonator formed by a mi¬ 
cropillar with moveable Bragg reflectors and a thin dielectric membrane at the 
node (or antinode) of the resonator. The mechanical displacement of the mem¬ 
brane quadratically couples to the cavity photon number. In addition, we have 
a monochromatic laser held with frequency ul applied to drive the cavity. The 
corresponding Hamiltonian for the scheme is given by, 


H = + (1) 

where Hq and H^i are given by Eqs. (2) and (3). In our model, H^i is the part 
where all the two subsystem cavity photons, phonons (due to mechanical motion 
of the membrane) respectively interact with their corresponding reservoirs which 
are collection of harmonic oscillators. We write 


Hq — LOc^^ a + + 9 opt 


a 


+I2[, 






( 2 ) 


where a and h correspond to held operators for cavity photons and phonons (due 
to mechanical motion of the membrane) with frequencies uoc and uom^ respectively. 
The third term in R.H.S of Eq. (2) describes the quadratic optomechanical coupling 
with strength gopt between the cavity held and the mechanical motion of the 
membrane and the fourth term describes driving process of the cavity held, with 
Q is the Rabi frequency of the driving held. Meanwhile, Hji takes the following 
form: 


for cavity mode 

Hr = 


)+s: 


for moving membrane 
<^k'n\,nki + '^gk' (^[,6 + 6 ^%'), 
k' k' 


(3) 

where the first two terms on the R.H.S of Eq.(3) represent the damping of cavity 
mode through reservoir of harmonic oscillator operators mj^ and while the third 
and fourth terms are responsible for damping of membrane through the harmonic 
oscillator operators nj^f and • In the rotating frame with the driving frequency 

(jOl^ the Hamiltonian Hq reduces to 


a^a I 


Hq — a + + gopt 

where Ac = {cjc 
quency ojl. Expanding the quadratic term 
Hamiltonian Hq in the following form: 


T 12 T 


(4) 


ul) is the corresponding detuning from the driving laser fre- 

21 


9opt 


a^a 


+ b 


')■ 


we have the 


Ho = Z\ca'*’a + + gopt [h h (P^ + P + 2pb + 1 


^ (5) 
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Coupling strength gopt should satisfy the condition {ujm + ^sgopt) > 0 for the sta¬ 
bility of the membrane m, where s is the number of photons inside the cavity. 
The corresponding Heisenberg equations of motion for the field operators of dif¬ 
ferent subsystem (cavity, membrane) are given by Eqs. (6) and (8), whereas the 
corresponding equations of motion for reservoir operators are given by Eqs. (7) 
and (9): 


d ^ ^ 

— a — —i/AcCi — igopt 


l^d (P^ + + 2pb + j — — i ^ dk'^k • 

k 


(6) 


d ^ ^ ^ 

—TU/j. = — igkd- 

^6 = -iwMb - 2igopt [a^o ( ^ + 

k' 

, - 


( 7 ) 

( 8 ) 

( 9 ) 


Now the corresponding operator equations for different reservoir need to be ex¬ 
plicitly solved further to remove the terms of reservoir operators in Eqs (6) and 
(8). Integrating Eq. (7), one gets 


mfc(i) = rnfe(0)e f dt d{t )e (10) 

Jo 

In Eq. (10), the first term represents the free evolution of the reservoir modes, 
whereas the second term arises from their interaction with the cavity. The reservoir 
modes 77ij^(t) can be eliminated from the Eq. (6) by substituting the value of m/c(t) 
from Eq. (10) in Eq. (6). So, we have 


= -iAcd-igopt (P^ + + 2pb + l)j -iH-y^ gl j dtd{t)e ^ ^+Ea(t), 

( 11 ) 

where Fa{t) = gk'f^ki.^)^ is the noise operator for the cavity which 

depends on the reservoir variables. The term containing the integral is expressed 
as 

/ di a(i ^ Taa(i). (12) 

fe 7o 2 

Thus, from Eqs. (11) and (12), we obtain 


^g = -iAcd - igopt [h (^P^ + P + 2pb + 1^j - iQ - ^Fa d{t) + Fa{t), (13) 

where Fa is the damping constant for cavity mode. Similarly, we have equation for 
mechanical motion of the membrane as follows: 

= -iiOMb - 2igopt [d^d (b + b + Fb{t)-, (14) 

where Ft,(t) and are noise operator and damping constant for moving membrane 
respectively, given by 

Fbit) = -*X]5fe'dfc'(0)e“*“'='‘, 
k' 


( 15 ) 
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and 






(16) 


2.1 Quantum Correlations 

We have obtained the coupled equations involving the mean of the operators 
a and b , their corresponding adjoints and their odd and even products. Since 
both the reservoirs are big as well as always in thermal equilibrium, we have 
{Fa{t)) = {Fb{t)) = 0. We also have, (i)S(i)) = (F^(i) d{t)) = (^Faitp(t)') = 
(t) a(t)) = 0. Similarly, quantum correlations between any subsystem (cavity, 
moving membrane) with different noise operator other than its own coupled reser¬ 
voir’s noise operator vanish altogether. The following set of equations has been 


obtained by using the Heisenberg-Langevin approach given in m- 

^ (a) = -iAc (a) - igopt ^ (^db^b^ + (a)] 

-iQ-^Fa(d). (17) 

^ = iAc + igopt + 2 + (a^)] 

+il2-(18) 

^ (b) = -icoM (b) - 2igopt [(a^aS) + _ Ipb (b) . (19) 

^ (b^'^ = iujM (b^^ + ‘2igopt [(a^ab^^ + {^a^aS^j - ipf, . (20) 

A ^a^'a) = -iQ ((a'!') - (a)) - Fa (Pa) + Fan^h- (21) 

P (Ph) = -2igopt [(a^a P'^) - (Pd p)'^ - Fb (Pb) + Fbph. (22) 

A (dP) =i{cuM- ^c) (ab^) -in(P)-^ (Fa + Fb) (dP) (23) 

-\-Wopt 1^2 + 2 ^^d^d^S^^ — — ^dS^^j . 

^ (^Pb^ = i (Ac - ujm) (^Pb^ + (b^ - ^ ( A + A) (^Pb^ (24) 


-\-i 90 pt + ^d^S^^ — 2 ^d^^dS^^ + 2 (^(^PPP"^ — ^d^^dS^^ + ^d^S^j . 

^ (db) = -i (Ac + ujm) (db) - iO (b) - ^ (Fa + Fb) (db) - igopt [s (dP) + (dP^b) 

+2 ^d^d^S^^ + ^dS^^ + 2 + ^d^d^S^^ + • (25) 
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^ = i (Ac + ujm) (^a'P'^ + iO (^P"^ - i (ra + A) (^PP"^ + igopt [2 + (^PPP"^ 

+2 + 2 + (^PP‘^b^ + ^d^^dS^^^ + . (26) 

P ^d^^ = -2iAc (^P'^-2iO {p-Fa (P^-2igopt + (PP^ + 2 ^d^S’^'S^ + (py^ • 

(27) 

^ (P‘^"^ = 2izAc ^d'*'^^+2zi7 ^d'*'^-ra (P‘^^^2igopt ^(^P^P‘^'^ + + 2 ^d'*'^S'*'S^ + • 

(28) 

P = -2iuM (P^ - A - 2z^opt |^{|d'*'d^ + 2 (^d'^'dS^S^ + 2 (^d^dS^^j . (29) 


P = 2zct;M (^P‘^'^+‘^Wopt |^<||d^d^ + 2 (^d'*'dS'*'S^ + 2 {||d'*'dS^^^j . (30) 

where are the thermal mean number of cavity photons and membrane 

oscillations respectively. 


2.2 Coupled Equations and Decorrelation of Higher-Order Operators 

The above set of equations are not closed, constitute higher order operator prod¬ 
ucts, and need to be solved numerically with approximations. We proceed to decor¬ 
relate all the higher-order correlations in the above equations. The above set of 
equations will then be closed up to second order when we apply the decorrelation 
method. We proceed to decorrelate the higher-(third-and fourth-) order quantum 
correlations present in the above equations like [28] , which studied the correlation 
of photon pairs from a double Raman Amplifier; a hybrid quantum optomechan¬ 
ical system containing a single semiconductor quantum well m- This approach 
corresponds to truncation of higher-order operator products in order to solve for 
all the second-order correlation functions. A similar kind of approximation has also 
been used in Ref. m to study the dynamics of a two-mode BEC beyond mean 
field approximation: 

(ABC) « [(A) {BC) + {AB) {C) + {AC) {B)] (31) 

and 

{ABCD) « [{AB} {CD) + {AC) {BD) + {AD) {BC)]. ( 32 ) 
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2.3 Temporal Evolutions 


After applying the decorrelation method, the above set of equations reduced to a 
closed set of coupled equations given in Appendix A. These equations contain all 
possible first-order operator averages and second-order correlation functions. The 
closed set of equations, given in Appendix A, is further solved numerically for the 
transient dynamics as well as nonclassical photon statistics. 

Figure 1 shows the temporal evolution of the mean number of excitations for the 
cavity photons and phonons due to mechanical motion of moving membrane 

for a set of different detuning parameters. For a detuning of the order of 
mechanical frequency cavity is driven around red sideband associated with 
the membrane oscillator and we have photonic oscillations (^d^d^ with phonon 

sidebands Both mean photonic and phonon excitations are asymmetric 

in shape also due to strong quadratic optomechanical coupling. As the detuning 
increases, cavity is driven far from red sideband regime due to which amplitude 
of the oscillations for are reduced as shown in Figure 1(c). For a very large 

detuned cavity, oscillations for is completely reduced to its initial value, 

whereas for we have a regular oscillations over the whole range of time scale. 

So, for a quadratic coupled optomechanical system cavity should be driven in red 
(blue) sideband regime of membrane oscillator. This is also required for generation 
of sub-Poissonian light from the driven cavity in our system Hamiltonian. 

In Figure 2, we observe the effect of cavity decay on the temporal dynamics 
of the system. In the weak cavity decay regime, we have undamped oscillations 
for photonic as well as phonon oscillations as shown in Figure 2(a). However, in 
the strong cavity decay regime, amplitude of oscillations for both as well 

as show damping with time as shown in the Figure 2(b). Although, in both 

the cases decay rate for phonons remains the same. For strong cavity decay, 
photons leak out of the cavity very rapidly and this leads to damped oscillations 
for mean phonon number also with time. 

We also study the temporal dynamics for finite thermal phonon numbers in 
Figure 3. We can see that the mean phonon number due to mechanical 

motion of membrane increases and saturates with time due to thermal noise (for 
finite n^^) even for the strong decay regime of cavity as well as mechanical mode 
as shown in Figure 3(b). So, a finite thermal noise due to phonons is found to 
increase the mean number of excitations for the dielectric membrane. This is due 
to dependence of mean excitations number of phonons on thermal noise 

and can be also seen in Eq.AG. In this case mean excitations number for photonic 
and phonons also oscillate in opposite phase upto certain range of time. 
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Fig. 1: (Color online) Mean number of cavity photons (blue line) and 

phonons of moving membrane (green line) for chosen set of parameters 

9optl^M = 1 - 4 ; Fa/LOM = 0 . 01 ; Fii/lom = 0 . 001 ; O/ujm = 0 . 6 ; = 0 ; 

where ujm is frequency of mechanical motion of the membrane chosen to 1 in our 
numerical simulations, a) AcjoJu — 0-5; b) Adoou = 1-0; c) AcjoJu — 2.0; d) 
— 5.0. 
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(a) (b) 

Fig. 2: (Color online) Mean number of cavity photons (blue line) and 

phonons of moving membrane (green line) for chosen set of parameters 

= 1-0; QoptloJM = 1-4; FijIum = 0.001; O/um = 0.4; = 0; a)For 

weak cavity decay Fajo^M — 0.01; b) For strong cavity decay FajojM = 0.1. 




Fig. 3: (Color online) Effects of thermal phonons on the mean number of cavity 
photons ^d^d^ (blue line) and phonons of moving membrane (h^b^ (green line) for 
chosen set of parameters AcjuJu — 1-0; Qoptl^M = 1-4; FljujM — 0.4 under strong 
decay regime FajoJu — FijIojm = 0.1; a) For = 0.0; b) = 2.0. 
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3 Photon Statistics 


In addition to the mean excitation numbers for cavity as well as mov¬ 
ing membrane we have computed the normalized two-boson correlation 

functions for the cavity field quadratically coupled membrane gl^\o) , and 

the cross-correlation between them using |27[I28[I29) . 


9?Ho) 


(dtd)^ 


2 ^d^d^ ^d^d^ + 

(dtd)^ 


(33) 



2^S‘^S^ (b^'^ + 

{my 



^d^S^ ^S^d^ + ^d^d^ + (bd^ 

(dtd) 


(34) 


(35) 


If g^\o), gl^^\o) and g^^^ (0) satisfy the inequality 5'^^(0) < 1, X = a,b,ab, 
then the statistics of the bosonic systems are referred to as sub-Poissonian. Statis¬ 
tics with ^^^(0) = 1 and ^'^^(O) > 1 are similarly referred to as Poissonian and 
super-Poissonian, respectively. Here, we have studied the temporal dynamics of 
second-order correlation functions, namely the self-correlation and cross correla¬ 
tion for cavity and oscillating membrane modes. 


In photon blockade effect, where coupling of a single photon to the system 
hinders the coupling of the subsequent photons, we have 5 'i^^( 0 ) < 1 for the cavity 
mode. Similarly, in photon induced tunneling regime, coupling of initial photons 
favors the coupling of the subsequent photons and leads to condition 5 'i^^( 0 ) > 1. 
We examine the effects of various parameters on the temporal behavior of 5'i^^(0), 
^^^^(0) and 5 'q 5 ^( 0 ) as well as under different decay regimes for cavity photons as 
well as phonons due to micromechanical motion. 


The second order correlation functions, namely the self- and cross-correlations 
for photonic and membrane oscillations are shown in Figure 4. We study the tem¬ 
poral dynamics of g^\o), gl^^\o) and for different cavity detunings Ac. 

It can be seen that for a resonant cavity i.e. Z\c = 0 and for a shorter period of 
time ~ 2, all the three correlations follow sub-Poissonian photon statistics as 
shown in Figure 4 (a). As UMt increases, ^i^^(O) tends to become super-Poissonian, 
whereas g^^^ (0) oscillates from sub-Poissonian to super-Poissonian. For a finite de¬ 
tuning Ac ~ ojm, g^\^) as well as 5 '^^^^( 0 ) both oscillate periodically from sub- 
Poissonian to super-Poissonian whereas g^^ (0) mostly remains sub-Poissonian as 
shown in Figure 4(b). For a very far off detuned cavity, g^\o) oscillates from Pois¬ 
sonian to super-Poissonian and it never becomes sub-Poissonian, whereas gl^^\o) 
behaves almost the same and oscillates from sub-Poissonian to Poissonian 
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as shown in Figure 4(c) and Figure 4(d). For a very large cavity detuning, a sin¬ 
gle photon can not be resonantly excited into the cavity, while the probability of 
finding two or more photons resonantly enhanced m- So, we do not have photon 
blockade effect for larger cavity detuning at any period of time UMt. 

Similarly, for a finite detuning Z\c, we study the effect of optomechanical cou¬ 
pling strengths gopt on various correlations in Figure 5. For the case gopt compa- 
rable to driving field i7, ga "^^0) and g^ "^^0) oscillate from sub-Poissonian to 
super-Poissonian, whereas g^^\o) remains mostly sub-Poissonian except at some 
points as shown in Figure 5(a) and Figure 5(b). For a very strong optomechanical 
coupling gopt^ oscillates from Poissonian to super-Poissonian, whereas 

remains same and varies to Poissonian limit as shown in Figure 

5(d). 

Furthermore, we examine the effects of photonic and phonon decay rates on 
various two-boson correlations in Figure 6. For very weak decay rates of cavity 
and membrane oscillation, ^i^^^(O) varies from Poissonian to super-Poissonian and 
never becomes sub-Poissonian. Although ^'^^^^(O) oscillates from sub-Poissonian to 
super-Poissonian and g^^\o) remains mostly sub-Poissonian as shown in Figure 
6(a). For a strong cavity decay rate Fa, becomes sub-Poissonian over a 

large scale of time and even becomes smaller than 0.5 as shown in Figure 6(b). 
This means, we have a photon blockade effect only in resolved sideband regime 
< 1 but not in the deep-resolved-sideband regime PalooM << 1 • While, 
for a very strong cavity as well as membrane oscillation decay rates, all the three 
correlations decay with time rapidly as shown in Figure 6(c). So, strong photonic 
and phonon decay rates play a positive role for getting sub-Poissonian photon 
statistics from a driven cavity in our system Hamiltonian. 



12 


Shailendra Kumar Singh, S. V. Muniandy 





(c) (d) 

Fig. 4: Second-order autocorrelations (0) (black solid line) for cavity mode, 
(0) for moving membrane (red dashed line) and g^^ (0) (blue dotted line) 
corresponding cross-correlation between photons and phonons for chosen set of 
parameters govij^M = 1-5; ra/uM = 0.01; = 0.001; = 0.6; = 

= 0; a) AcjoJu = 0.0; b) Acjuju = 1-3; c) Ad^u = 2.5; d) AcjuJu = 4.0. 
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(a) (b) 




Fig. 5: Second-order autocorrelations (0) (black solid line) for cavity mode, 
(0) for moving membrane (red dashed line) and g^^ (0) (blue dotted line) 
corresponding cross-correlation between photons and phonons for chosen set of 
parameters Acjuju = 0-5; Fajoju — 0.01; = 0.001; Fljuju — 0.6; = 

= 0; a) govtl^M = 0.8; b) govtl^M = 1-7; c) govtl^M = 3.0; d) govtl^M = 6.0. 
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(a) (b) 



(c) 


Fig. 6: Second-order autocorrelations (0) (black solid line) for cavity mode, 
(0) for moving membrane (red dashed line) and g^^^ (0) (blue dotted line) 
corresponding cross-correlation between photons and phonons for chosen set of 
parameters AcjoJu — 1-0; Qoptl^M = 1-4; ^Iojm = 0-4; = 0 under 

different decay regimes; a) FajojM = 0.01, — 0.001; b) FajuM = 0.1; 

^hl^M — 0.001; c) FajcOM = ^hl^M — 0.1. 
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4 Conclusion 

We have studied a qudratically coupled optomechanical system formed by a mi¬ 
cropillar with moveable Bragg reflectors and a thin dielectric membrane at the node 
(or antinode) of the resonator through the Heisenberg- Langevin approach. We 
have obtained coupled equations up to the second order without any approxima¬ 
tion. We decorrelate the higher-order correlation functions in the coupled equations 
to obtain a closed set of equations. This enables the study of temporal dynamics 
of the optomechanical system as well as two-boson correlation functions under dif¬ 
ferent scenarios. In the resolved sideband regime, mean number of cavity as well as 
mechanical excitation show oscillations. Furthermore, amplitude of these oscilla¬ 
tions changes as we change the cavity detuning. We have also studied the temporal 
dynamics under the different decay regimes of cavity and membrane oscillations. 
We have discussed the effect of thermal noise due to phonons on the dynamics 
of the optomechanical system. Furthermore, we analyzed the temporal dynam¬ 
ics of the second-order correlation functions, namely the self and cross-correlation 
of the cavity and mechanical modes under different cavity detunings as well as 
optomechanical coupling strengths. We obtained the regimes where all the three 
correlations display strong sub-Poissonian photon statistics particularly at very 
small cavity detuning as well as strong cavity decay rate. Our study is useful for 
coherent control of photon statistics as well as photon and phonon correlations in 
quadratically coupled optomechanical systems. 
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A Decorrelated and Closed set of Coupled Equations 

^ (a) = -lAc (a) -iQ - ipa {a) - igopt + 2 (^^)} 


(36) 
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^ = iAc + ii? - + igopt [{(^^) + 2 {^0} 

+ {(at) (S 2 ) + 2 (ats) (S)} + 2 {(at) (StS) + (atSt) (S) + (atS) (St) j + (at)]. 

^ (s) = -i(^M (S) - (S) - 2igopt [(at) {(aS) + («^^)} 

+ (ata) {(S) + (St)I + (a) {(«^^) + 

^ (St) = iujM (st) - (st) + 2igopt [(«^) {(^^) + {^^0} 

+ (ata) {(S) + (St)| + (a) {(atS) + • ( 39 ) 

^ (ata) = -in ((at) - {a}) - Fa (ata) + Fan^^. ( 4 o) 

^ (StS) = - 2 igopt [{(ata) (St^) + 2 (atSt) (aSt)| - {(ata) (S^) + 2 (atS) («^)|] 

-n(StS)+r,n,\. ( 41 ) 

^ (aSt) = i {(jjm - Ac) (aSt) - in (St) - 1 (Ta + r^) (aSt) 

+igopt [2 {2 (ata) (aS) + (atS) (a^) j - {(«^^) (^^) + 2 (aS) (^^^)} 

+2 {(2 (ata) (aSt) + (atSt) (a^)) - (2 (StS) (aSt) + (aS) (^^^))| 

-(aSt)( 3 (St 2 ) + l)]. ( 42 ) 

^ (atS) = i {Ac - (jjm) (atS) + in (S) - i (Ta + a) (atS) 

+wopt [{(^^^) ^ (atSt) (StS)I - 2 {(at^) (aSt) + 2 (ata) (atSt) j 

+2 {((atSt) (S^) + 2 (atS) (StS)) - ((at^) (aS) + 2 (atS) (ata))} 

+ (atS) (3 (^^) T t-)j. (43) 


— (aS) = —i {Ac + (^m) (aS) — zi? (s) — {Fa + -T^) (aS) 

dt ' ' '/ '/2 '/ 

[2 (aSt) + 3 (aS) ((S^) +1) j - [(^^) ^ (aSt) (StS)j 

-2igopt [(2 (ata) (aSt) + (atSt) (a^)) + ((«^^) (^^) + 2 (aS) (^^^)) 

+ (2 (ata) (aS) + (atS) (a^))j. (44) 

^ (atSt) = i{Ac fujm) (atSt) + zi 7 (St) - i (Ta TA) (atSt) 

Figopt [2 (atS) + 3 (atSt) ((St^) + i)j + igopt [(^^^^) + 2 (atS) (StS)j 

+2igopt [((^^^) + 2 (ata) («^^)) + («^^) + 2 (atSt) (StS)) 

+ (2 (ata) (atSt) + (aSt) (at^))j. (45) 
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— (o^) = —2iAc (o^) - 2iQ (a) — Fa (a^) - 2igopt (a^) 


(46) 


‘^Wopt 


(o,^) + 2 + 2 +2(db)(db 


^ = 2iAc ^a^^ - Fa + ‘^Wopt 


H“ ‘2‘igopt 


(at2) ((6t2^ + ^^2) + 2 (fets)) + 2 ( (atfet)" + (atfe)" + 2 (a 


— (h^^ = -2i(jjM (b^^ - Fb - 2igopt ^a^a^ (48) 

-^igopt l^^a^a^ ((^^) ^ ^ ^ ^ ^ • 


— = 2zcaM “ A + 2i5fopt {^a^a^ (49) 

Fdigopt l^^a^a^ ^ ^ ((^^) + ^ ^ • 



